Effects of edge disorder in nano-scale antiferromagnetic clusters 
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We study the distribution of local magnetic susceptibilities in the two-dimensional antiferromag- 
netic S = 1/2 Heisenberg model on various random clusters, in order to determine whether effects 
of edge disorder could be detected in NMR experiments (through the line shape, as given by the 
distribution of local Knight shifts) . Although the effects depend strongly on the nature of the edge 
and the cluster size, our results indicate that line widths broader than the average shift should be 
expected even in clusters as large as ~ 1000 lattice spacing in diameter. Experimental investigations 
of the NMR line width should give insights into the magnetic structure of the edges. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Cx, 75.40.Mg 
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Experimental, theoretical, and computational studies 
of two dimensional (2D) quantum spin systems have 
taken a prominent place in condensed matter physics 
during the past two decades. The 5 = 1/2 Heisenberg 
model has successfully explained the bulk antiferromag- 
netism in the parent compounds of the high-T c cuprates 
P, 0], and other quasi-2D materials 0, It na s also 
been realized that controlled studies of systems with im- 
purities can give additional valuable information on the 
electronic structure and interactions in these strongly- 
correlated materials Impurity effects are also of in- 
terest in their own right, as they reflect fascinating quan- 
tum phenomena not present in translationally invariant 
systems 

In this Letter, we address the physics of defects be- 
yond single-spin impurities in 2D antiferromagnets, by 
considering various forms of edge disorder. Prominent 
effects of effectively free chain ends are known in ID 
systems, where detailed comparisons of theories, com- 
putational model studies, and experiments are possible 
[lfj,lll|. In contrast, 2D edges have not been paid much 
attention to, presumably because one would naively ex- 
pect their influence to be relatively small in experiments 
probing spatially averaged properties, due to typically 
small edge-to-bulk ratios. With the increasing focus on 
physics on the nano-scale, it should also be interesting to 
investigate small antiferromagnetic clusters, where edge 
effects could dominate the physics. To this end, an ini- 
tial quantum Monte Carlo (QMC) study of both smooth 
and rough edges in the 2D 5 = 1/2 Heisenberg model 
was recently carried out [l2| ■ A suppression of the mag- 
netic susceptibility at a smooth edge was found at low 
temperature (T), contrary to the naively expected en- 
hancement due to the smaller coordination number of 
the edge spins. The edge contribution to the susceptibil- 
ity is logarithmically divergent for T — ¥ 0, as was later 
found also within a continuum field-theory description 
[HI . Another interesting observation was a dimerization 



pattern at the edge, which can be seen as a precursor 
to valence-bond solid state that can be realized if addi- 
tional interactions are included 14 1. These effects may 



not be easy to observe experimentally, however, because 
edge roughness masks this behavior [12|. The roughness 
introduces subsets of spins which are effectively weakly 
coupled to the bulk system, leading to a strongly en- 
hanced susceptibility. The behavior reflects a complex, 
and still not well understood, interplay between geomet- 
ric roughness at the microscopic scale and the collective, 
macroscopic behavior of interacting spins. 

The nature of the clusters in powder samples of quasi- 
2D antiferromagnets such as La2CuC>4 is currently not 
known precisely, and, to our knowledge, there have not 
been any efforts to investigate possible effects of finite- 
size clusters in, e.g., nuclear-magnetic-resonance (NMR) 
experiments. Cluster diameters > 1000 lattice spacings 
can be expected [lj|- One would expect free edges to 
lead to broadening of the NMR line (a distribution of 
Knight shifts), as in ID jlo| . The question is, whether 
this broadening can be observed, and what information 
it can provide on the structure of the edges. 

We here present a systematic study of the Knight-shift 
distribution in finite 5 = 1/2 Heisenberg 2D clusters con- 
structed with varying amounts of edge roughness. We 
calculate the width of the distribution of local magnetic 
susceptibilities as a function of the cluster size. The re- 
sults indicate that edge effects should be very prominent 
even for relatively large clusters, of width 100-1000 lattice 
spacings, at temperatures where the cuprates are para- 
magnetic (i.e., above the Neel temperature, Tn, where 
order sets in due to weak 3D couplings or anisotropics). 

With only on-site hyperfine couplings taken into ac- 
count, the NMR Knight shift of a Cu nuclear spin at 
lattice site r is proportional to the local susceptibility 
Xi(r), which is given by 



Xi( r ) 
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where /3 — 1/T (in units where kg = 1). In an infinite 
uniform system xi( r ) = Xi t ne bulk susceptibility. Dis- 
order or open edges in a finite system lead to local vari- 
ations and, thus, a broadened NMR line. In reality, in 
cuprates there are significant transferred hyperfine cou- 
plings also to nearest-neighbor Cu sites and the Knight 
shift should be modified accordingly [l6|. Here, in this 
initial proof-of-concept study, we do not consider these 
couplings and just report results of the completely local 
susceptibility. The methods we use can, however, easily 
be extended to any NMR form factor. 

Distributions of susceptibilities due to isolated vacan- 
cies (corresponding to Cu substituted by Zn in cuprates) 
in the 2D Heisenberg model 17, 3| and ladder systems 
[HI have previously been studied using QMC simula- 
tions. Here we use similar techniques to study open- 
boundary clusters with various types of edges. We use 
the stochastic series expansion method with efficient loop 
updates and improved estimators flil . [2fj| to evaluate ([Tj) 
for various types of clusters at different temperatures. 

We build three different types of clusters: A) As in 
Ref. 12[, starting from an open L x L lattice, we tra- 
verse the 4(L — 1) boundary sites and remove each spin 
with probability p — 1/3 or couple a new spin to it with 
the same p, doing nothing with probability 1 — 2p. B) 
With spins randomly occupying the sites of an open Lx L 
lattice with probability p, we identify the largest cluster. 
C) Starting from a single occupied site on an infinite lat- 
tice, we add neighbors to it with probability p, filling the 
neighbors of added spins with this same probability, until 
a cluster of a desired size N has formed. We here only 
consider the case p = 0.7 and 0.6 for B and C clusters, 
respectively. Since we want to isolate the effects of edges, 
in all cases we also fill any internal vacancies, so that all 
sites within a single boundary are magnetic. The C clus- 
ters are close to the percolation point (p w 0.59) and are 
therefore very (unrealistically) irregular in shape, thus 
serving as an extreme case. 

In all cases, we here initially consider only clusters 
with the same number of spins on both sublattices. This 
is done in order to avoid a trivial enhancement of the 
susceptibility originating just from the fact that clusters 
with sublattice imbalance have ground states with non- 
zero total spin. On the other hand, in real systems, sub- 
lattice imbalanced clusters will of course be present, and 
internal vacancies or defects can be expected as well. Our 
calculations here are intended to better isolate the effects 
solely due to the edges. We will also investigate the ef- 
fects of the fluctuation of the total ground-state spin in 
the latter part of this paper. 

Fig. Q] shows the local susceptibility landscape of an in- 
tact open lattice as well as a rough-edge cluster of type A. 
For the smooth-edge cluster, the corners have the largest 
response. Because of this, in the presence of a weak ex- 
ternal field the neighbors of the corners would experience 
an effectively negative field, and, thus, the susceptibility 




FIG. 1: (Color online) Local magnetic response for two sys- 
tems; an intact 32 x 32 lattice (left) and a rough-edge A clus- 
ter (right). The upper and lower graphs are for (5 = 2 and 
/3 = 10, respectively. The intensity color-codings correspond 
to lowest-highest response according to: 0.056 < \i < 0.186 
(left, upper), -0.19 < Xi < 0.47 (right, upper), -0.087 < 
Xi < 0.177 (left, lower), -1.8 < xi < 2-9 (right, lower). 

of these spins is negative at low temperatures. This kind 
of staggered susceptibility pattern continues away from 
the corners on a length-scale which should be related to 
the exponentially divergent correlation length Q ■ A sim- 
ilar effect around an isolated impurity has been studied 
in great detail previously [17| • Note that the low-T sup- 
pression of the edge susceptibility found in [ij} refers to 
the uniform component (average over spins at a fixed dis- 
tance from the edge), whereas the staggered response is 
enhanced at the edge. 

Our main interest here is in the response at and close to 
a rough edge, where strong staggered susceptibility pat- 
terns can appear around some spins for the same reasons 
as discussed above. Because of interference among sev- 
eral randomly located high-response centers, these effects 
are amplified in some regions and damped in others, as 
seen in the type A cluster in Fig. [TJ leading to a complex 
behavior with much larger fluctuations in local suscepti- 
bilities than in the intact open system. 

To characterize the susceptibility distribution as a 
function of the cluster size, we define the average radius 
of a cluster c with respect to the site closest to the "grav- 
itational center" rg, and further average this radius over 
different cluster realizations; 



(2) 



Here r° is the position of an edge spin and is the 
number of edge spins in cluster c. 

Fig. [2] shows histograms of the local susceptibility (cor- 
responding to the NMR profile) at inverse temperatures 
(3 = 2 and 4. Note that, in the context of cuprates, 
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FIG. 2: (Color online) Local susceptibility distribution for 
clusters of types A, B, C, as well as intact open lattices, 
all at two different inverse temperatures; /3 — 2 and /3 — 4. 
Note that different scales are used for the intact lattices (of 
sizes L = 32, 128, 512 corresponding to the average radii R 
indicated). The results for A, B, and C clusters are averages 
over 200-1000 different random realizations. 

these temperatures are rather high, with ft = 4 corre- 
sponding approximately to room temperature (close to 
Tn for La2CuC>4). NMR experiments have, however, in 
the past been carried out even up to T = 900 K [2l| . 
and the behavior there is still in good agreement with 
the 2D Heisenberg model (22)]. For all clusters, the shift 
distribution is much wider at ft = 4 than at ft = 2. At 
the higher temperature ((3 — 2), the bulk susceptibility 
rapidly starts to dominate as the cluster size is increased. 
At j3 — 4 there are, however, still very significant broad- 
ening effects for clusters as large as R ~ 300. Experi- 
mentally, for fine-powdered cuprate samples, one would 
expect clusters roughly of this size or somewhat larger 
[15l | . For small intact L x L clusters at (3 = 2, the edges 
give rise to a separate large peak at xi ~ 0-1- There is 
also a much smaller peak due to the corners at higher xi 
(not seen in the figure). At j3 — 4 several smaller peaks 
originate from frames of spins at and close to the edge 
[l2j. These effects persist for larger clusters but become 
difficult to discern because of their low relative weight. 
Overall, the random edges lead to much broader peaks. 

To quantify the width of the susceptibility distribution, 
we calculate the standard deviation with respect to the 
average local susceptibility over all clusters, 
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^^rjEEW-W) 2 . (3) 
where c labels the individual clusters, N c is the number 
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FIG. 3: (Color online) Standard deviation of the local suscep- 
tibility distribution for the three types of clusters as functions 
of the effective radius (R) at several inverse temperatures. 

of occupied sites of the clusters, and r is a site label for 
these magnetic sites; thus, Xi( r ) ^ s the local magnetic 
susceptibility of site r on cluster c. Its average over all 
spins of all clusters is denoted (xi)- 

Fig. [3]shows results as a function of the average cluster 
size (R) for inverse temperatures in the range /3 € [1, 10]. 
In all cases, the line width first grows with the cluster 
size and then decreases. For any finite /3, as (R) — > 00 
we must have a x — ¥ 0, as the edge-to-bulk ratio vanishes 
for an infinite cluster and the edge effects can only extend 
a finite distance away from the edge at any T > 0. The 
initial increases in a x reflect the tails of the distributions, 
which only develop fully for large clusters. The tails are 
very significant at low temperatures, i.e., the rough edges 
influence the response far inside the bulk of the clusters 
(reflecting the exponentially divergent correlation length 
for T ~> 2]). At the lowest temperature {fi = 10) the 
line widths still increase for the largest clusters we have 
studied. Comparing the line-widths of the three types of 
clusters, we see that type C always produces the broadest 
lines (for given temperature and cluster radius). These 
cluster are also the ones that are the geometrically most 
rough ones, being constructed close to the geometrical 
percolation point. Real cuprate clusters are likely not as 
rough. Note, however, that even clusters of type A, for 
which the geometrical edge disorder is very shallow, lead 
to distributions only a factor of ~ 2 narrower. 

For all the disordered clusters with (R) « 300, the line 
widths exceed the average Knight shift for /3 = 3 ~ 5; 
relevant for cuprates in their paramagnetic state. At 
P = 4, A clusters with (R) = 295 give a x = 0.103, while 
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FIG. 4: (Color online) Size dependence of the line width a x 
for type A clusters constructed under the condition of equal 
numbers of spins on the two sublattices (filled symbols), and 
without this restriction (open symbols). 

the average shift is ( X i) = 0.0576, type B clusters with 
(R) = 294 give a x = 0.105 and ( X i) = 0.0571, and type C 
clusters with (R) = 301 have a x = 0.179, (xi) = 0.0622. 
Line broadening at this level should be clearly visible ex- 
perimentally (although the long tails of the distributions 
may be partially drowned by experimental noise). 

We now discuss the case of averages over clusters with- 
out the restriction of equal numbers of spins on both 
sublattices. Fig. [4] shows the line widths for the two en- 
sembles of type A clusters. At high temperatures, the 
line widths are almost equal (indistinguishable in the 
figure at /3 = 2) widths. At lower temperatures, the 
unrestricted ensemble gives significantly broader lines, 
however (about a factor of two at (3 = 6). For very 
large clusters we expect no differences, because the rel- 
ative fluctuation of the difference in sublattice occupa- 
tion is ~ 1/ \/W c . The convergence of the two ensembles 
with increasing (R) is seen clearly for j3 = 4. Based on 
these results, we expect the contributions from the sub- 
lattice imbalance to be small relative to the dominant 
boundary-roughness effects in the size and temperature 
regimes most relevant to cuprates. 

In summary, we have shown that significant NMR line- 
broadening should be expected due to rough edges of 
nano-scale antiferromagnetic clusters. In the tempera- 
ture range of relevance, e.g., for the cuprates in their 
paramagnetic state, the line widths produced by all the 
cluster types studied are larger than the average Knight 
shift, even for clusters several hundred lattice constants 
wide. The effects should therefore be experimentally ob- 
servable in fine-powdered samples. 

In light of these results, it may seem surprising that 
no anomalous broadening has been noted in experiments 
so far (perhaps suggesting that the edges actually are 
rather smooth) . It should therefore be interesting to sys- 
tematically search for edge effects in NMR and other ex- 
periments. In order to make the NMR broadening more 
prominent and investigate the size dependence, it would 
then be desirable to prepare powder samples with very 
small clusters, down to m 100 lattice spacings across. 



Current techniques (which are normally not intended 
to produce extremely small clusters) likely give clusters 
roughly an order of magnitude larger [15j . but other 
methods could perhaps reach smaller sizes. Systemati- 
cal studies of antiferromagnetic correlations and response 
functions in nano-scale systems may give further insights 
into the microscopic interactions and collective quantum 
phenomena in strongly-correlated electron systems. 
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